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TWO FACTORS INFLUENCING TEMPERATURE DISTRIBUTIONS 
AND THERMAL STRESSES IN STRUCTURES 


By William. A. Brooks, Jr., George E. Griffith, 
and H. Kurt Strass 


SIMMARY 


The influence of Joint conductivity and internal radiation on tem- 
perature distribution and thermal stresses has been discussed. Joints 
of poor conductivity can occur in normal fabrication procedure and 
greatly alter temperature distributions and increase thermal stresses. 
On the other hand, internal radiation tends to make the temperature 
distributions more uniform and thereby relieves thermal stress. 


INTRODUCTION 


Thermal stresses are unquestionably an important consideration in 
the design of supersonic aircraft. Such stresses are usually produced 
by nonuniform temperature distributions - the greater the temperature 
variation within the structure, the larger the thermal stresses. The 
present paper deals briefly with two factors which may affect the tem- 
perature distribution, and thus the thermal stresses, within a structure; 
these two factors are Joint conductivity and internal radiation. In 
order to indicate seme of the effects of these two factors without 
unnecessary structural complications, the basic structure considered 
consists of a length of skin with an integral or an attached web. In 
addition to theoretical results, experimental data were obtained by 
heating the structure either aero dynamically or by laboratory radiant- 
heat sources. 


SYMBOLS 


a coefficient of thermal expansion 

c specific heat 


E 


modulus of elasticity 
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e 

emissivity 


* 

h 

heat-transfer coefficient 


k 

h j 

joint conductivity coefficient 


— 

k 

thermal conductivity 



M 

Mach number 


— 

a 

heating rate 



t 

skin thickness 



T 

temperature 


w 

ib 

initial temperature 

-■ 


*1 

maximum skin temperature 



T 2 

minimum web temperature 


- 

■*aw 

equilibrium temperature 

- 


T 

time 


V 

w 

specific weight 



a 

thermal stress 



Subscripts : 



RC 

radiation and conduction 



C 

conduction 



MAX 

maximum 




DISCUSSION 




Simple Integral Structure 


- 

Consider first the simplest possible type of skin-web combination, 
one without joints and without internal radiation effects. Figure 1 
shows such a structure, a symmetrical integral H-section with a flange 
(or skin) of thickness t symmetrically heated at a constant rate qj 

* 
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within this structure heat transfer takes place by conduction only. The 
actual temperature T minus the initial temperature To, multiplied by 
the material thermal conductivity k and divided by the heat in g rate q 
and skin thickness t, is plotted against a time parameter in which t 
is the elapsed time, c the specific heat, and w the specific weight. 
The choice of these d im ensionless parameters makes it possible to present 
on the same plot the results for a given cross section of different 
materials subjected to different heating rates. The calculated results, 
shown by the solid lines, are obtained by a simple analysis assuming 
constant material thermal properties and no heat losses. Point 1, which 
is farthest from the heat sink afforded by the web, will be the hottest 
point in the structure, whereas point 2, which is farthest fraa the heat 
source, will be the coolest. As a consequence of assuming no heat losses 
and a constant heat input, the average temperature of the structure is 
proportional to the mass of the section and varies linearly with time. 

In order to check the theory, three tests were made involving constant 
heating rates of 5* l6, and 4l Bfcu/ft^-sec. Each test is represented 
by a different symbol, and the experimental temperatures are plotted for 
various times until the maximum skin temperature reached 450° F. The 
parameters are such that the data for the highest heating rate, shown by 
the diamonds, appear at the extreme lower end of the curves, whereas the 
data for the lowest heating rate extend the full length of the curves. 

The test results are in fairly good agreement with the theory $ any dis- 
crepancy is due largely to the assumption of constant thermal properties 
used in the theory. 

At any point in the structure, the thermal stress is proportional 
to the difference between the average temperature and the temperature 
of the point in question. When multiplied by the appropriate material 
properties, these differences can be converted into thermal stresses, 
such as shown in figure 2. 

In this figure thermal-stress distributions are shown for both the 
low and high heating rates of 5 and 4l Btu/ft^-sec for the times when 
the maximum skin temperature reached 450° F. Hie skin stresses appear 
in region A and the web stresses, in region B. The solid lines represent 
thermal stresses calculated from measured temperatures shown in figure 1, 
and the symbols represent experimental stresses obtained from strain-gage 
readings corrected for differential expansion between the gage and speci- 
men. The agreement between theory and experiment is quite good. As 
would be expected, the higher heating rate allows less time for heat to 
be conducted into the web, thus provides larger temperature differences 
in the structure, and consequently yields larger thermal stresses. For 
the heating rate of 4l Btu/ft^-sec the maximum skin stress, shown at the 
extreme right, and the maximum web stress, at the extreme left, are more 
than twice the values shown for the rate of 5 Btu/ft^-sec, despite the 
fact that the skin temperature rise was the same. 
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Another significant point is that the maximum stresses are not 
related in the same proportion as the heating rates, which have an 8 to 
1 ratio. The stresses shown here are those which existed when the skin 
temperature T-^ was 450° P. If, however, the skin temperature were not 

limited to 450° F and a comparison were made at the same value of time, 
then, indeed, the stresses would "be proportional to the heating rate. 

It is also possible to present the stresses in a dimensionless form 
as was done for the temperatures . Figure 3 shows the dimensionless 
thermal stresses at various times during the progress of the three tests 
at heating rates of 5* l6, and 4l Btu/ft^-sec. It is again evident that 
the high-heating-rate data are grouped at the lower end of the theoretical 
curves which are represented by the solid lines, whereas the data for the 
low heating rate extend the full length of the curves. The theoretical 
stresses are obtained by using the theoretical temperatures of figure 1. 
The agreement between theory and experiment is about the same as that 
in figure 1, and any discrepancy is due to the seme reason. 

The results of these three figures indicate that fairly accurate 
temperatures and stresses can be predicted in a simple integral structure 
where only conductive heat transfer need be considered and, further, that 
strain-gage readings, if treated properly, can give reasonable stress 
results. However, if this structure had contained any joints, greater 
temperature differences would have resulted, and accurate calculations 
would have become somewhat more difficult. 


Structure With Attached Webs 

In order to show some of the effects of joints, the temperatures in 
one skin-web combination of a multiweb wing, a more realistic aircraft- 
type structure which contains riveted joints, will be examined. When 
six 20-inch-chord, multiweb wings identical in cross section - that is, 
with the same overall size, skin thickness, material, and so forth - were 
tested at sea-level conditions at a Mach number of 2 in a blowdown jet, 
the results were as shown in figure 4. This figure shows the temperature 
difference T-[_ - Tq between the maximum temperature in the skin (point 1) 

and the minimum web temperature (point 2) divided by the maximum possible 
skin temperature rise T av - Tq plotted against time. The upper dot- 

dash curve gives the theoretical temperature-ratio drop if no heat is 
conducted into the web, that is, if the joint conductivity hj is equal 
to zero, and the lower dot-dash curve gives the theoretical results 
for perfect joint conductivity. The experimental data for five of the 
wings, indicated by the shaded area, lie close to the curve for perfect 
joint conductivity. On the other hand, the data for the sixth wing show 
that the increased temperature drop is scmewhat closer to the curve for 
no joint conductivity. Since the temperature distribution depends mainly 



NACA TN 4052 


5 


upon heat conduction within the structure, this increased temperature 
difference can he attributed almost solely to a lower joint conductivity. 
The much poorer thermal conductivity of the joints of the one wing, can- 
pared with the other five wings, was not obtained by design but rather 
results from the small variations in otherwise identical structures 
expected fran normal fabrication techniques. The 50-percent increase 
in temperature drop of the one wing over the average of the other five 
wings also resulted in a substantial increase in thermal stress. In 
fact, the wing with the largest temperature drop fluttered and then 
suffered a dynamic failure; whereas, the one wing among the other five 
which was identical in every detail survived its test without damage. 
Construction details, such as the inclusion of ribs, invalidate further 
comparison with the other models. However, it seems obvious that joint 
conductivity can have an important, and sanetimes adverse, effect on the 
temperature variation in a structure and hence on the thermal stresses 
and, as a result, on the structural integrity. 

In figure 5 the maximum nondimens ional temperature drops have been 
plotted against a dimensionless joint conductivity parameter hjt/k 

where hj is the thermal conductivity of the joint. Calculated results 
are shown by the solid lines. In the calculations the value of 0.013 for 
ht/k, the Biot number (an index of the rate of external to internal, heat 
transfer), was determined by the experimental aerodynamic heat-transfer 
coefficient h. The symbols indicate data obtained from the tests and 
show a 6 to 1 variation in hj which apparently can be expected even 
from good shop practice. Data obtained by Barzelay, Tong, and Holloway 
(refs. 1 and 2) cover an even wider range in hj. Although a consider- 
able scatter in hj is shown by the five lowest test points, these 
points are in good agreement with the theory since they lie in a region 
where a fairly large change in joint conductivity has little effect on 
the temperature drop. Thus, the band for the five wings shown in fig- 
ure 4 masks a considerable variation In hj. The sixth point, which 
does not show such good agreement, is in a region where a small change 
in hj can alter the temperature drop appreciably. The spread of hj 
within the band shown in figure 4 is of the same order of magnitude as 
the variation from the top of the band to the curve for the one wing. 

If the maximum skin or web stress were plotted Instead of the tem- 
perature drop, a similar variation with hjt/k could be observed - the 

lower the value of hj, the higher the thermal stress, as evidenced by 
the failure of the one wing. 

Actual flight-test results are shown in figure 6 for a duplicate of 
one of the five wings which had very good joint conductivity. The Mach 
number and altitude for the flight are shown at the top of the figure. 

The temperature difference T-j_ - T£ between the skin and the web 
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center line is plotted against time. When this wing was attached to a 
rocket model, the results were as shown by the solid line, which lies 
much closer to the calculated curve for h^ = 300 Bfcu/ft^-hr-°F than 
to the curve for perfect joint conductivity. Since a value of hj = 300 
would lie slightly to the left of the test point for the one wing with 
the largest temperature difference shown in figure 5 > the joints of this 
wing would he considered poor frcm the standpoint of thermal conductivity. 
In this case failure of the wing was prevented by the inclusion of a 
chordwise rib. 


Internal Radiation in an Integral Structure 

As was mentioned previously, the effect of joints is to increase 
the magnitude of the thermal stresses. However, if the temperatures 
involved are sufficiently high, it is possible that internal radiation 
in a structure may bring about seme thermal-stress relief. The hotter 
parts of the structure transfer heat to the cooler parts by radiation 
and thereby lessen the temperature differences. Although there is an 
increased interest in this phenomenon, very little information can be 
found on this subject. 

Therefore, in order to evaluate some of the effects of internal 
radiation on the temperature distribution within a multicell beam, such 
as shown in figure 7, an analog solution was made to determine the 
temperature distribution of an infinitely long box whose internal surfaces 
have an emissivity of 1. The procedure was to assume a constant rate of 
heat input until the maximum skin temperature reached a prescribed limit 
and then to maintain this skin temperature while the web temperature 
approached the skin temperature. 

Heat was put into both skins at the rate of 5 Btu/ft^-sec until 
the max imum skin temperature T^ reached 1,200° F, at which time the 
skin temperature was kept constant. The maximum skin temperature Ti 
and the minimum web temperature Tg are plotted against time in seconds. 

The dashed lines are the solution when internal radiation is neglected; 
the solid lines represent the case when radiation is included. In this 
situation it can be seen that, when radiation is included, steady state 
is reached much more rapidly than when conduction only is considered. It 
can also readily be seen that the maximum temperature difference between 
the skin and the web, when conduction and radiation are considered, occurs 
at approximately 75 seconds and is about two-thirds of that when conduc- 
tion only is considered. 

In order to make a preliminary survey, other heating rates and skin- 
temperature levels were also investigated as shown in figure 8. Here, 
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in order to show one of the effects of internal radiation, the ratio of 
the maximum difference between skin temperature and web temperature for 
both radiation and conduction to the maximum, difference for conduction 
only is plotted against the skin temperature T-|_. For a given heating 

rate, the ratio of the temperature differences will approach a minimum 
value as the skin temperature increases. In the case of 1=1 Btu/ft 2 -sec, 
the limit is nearly reached at a skin temperature of 1,200° F. The curve 
for the highest heating rate will eventually reach the lowest limit but at 
a temperature far in excess of 1,200° F. For the conditions in the figure, 
radiation effects are most significant when heating rates are low and tem- 
peratures are high. 

Figure 9 shows seme preliminary experimental results . The specimen 

1 /I 

was made of — ^-inch Inconel with the cross section as shown (2 — inch- 

16 V 2 

square cells j and was 12 inches long. The emissivity of the internal 
surfaces was 7 about 0.8. The heating rates were approximately 5 and 
2-i Btu/ft 2 -sec. Although there is a measurable reduction in the maximum 

temperature difference, the reduction is not as large as indicated by the 
idealized theoretical results. 


CONCLUDING REMARKS 


In summary, two factors, joint conductivity and internal radiation, 
which influence temperature distribution and therefore thermal stresses, 
have been discussed. Even normal fabrication techniques can produce 
joints of such poor conductivity as to cause the temperature differences 
to increase markedly over those of an integral structure. On the other 
hand, at the higher temperatures internal radiation has the effect of 
making the temperature differences less severe and the thermal stresses 
smaller. At the higher temperatures these two effects tend to cancel 
one another. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 

Langley Field, Va., March 6, 1957- 
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TEMPERATURES IN RADIANTLY HEATED H- SECTION 
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Figure 1 


THERMAL STRESSES IN RADIANTLY HEATED H-SECTION 



Figure 2 
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Figure 5 


TEMPERATURE DIFFERENCE FOR MULTIWEB WING IN FLIGHT 
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Figure 6 
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TEMPERATURE HISTORIES FOR BOX BEAM WITH INTERNAL RADIATION 
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Figure 7 


INTERNAL RADIATION EFFECTS FOR VARIOUS HEATING RATES 



Figure 8 
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EXPERIMENTAL RESULTS FOR INTERNAL RADIATION 
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